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■ It is well established that intermediate (2-8 Mq) and high-mass 

(> 8 M©) stars form in cluster environments (e.g., Kurtz et al. Concerning the properties of the intermediate/high-mass 

2000; Evans et al. 2009), and that at the very first stages of their YSOs forming within the cloud, a few studies toward massive 

formation, they are embedded within dense and cold gas and star-forming regions have revealed that two objects formed in 

dust from their original natal cloud. However, it is not clear to the same environment may have very different properties in the 

what extent the initial conditions of the natal cloud are deter- ejection of matter, with one YSO driving a highly colhmated 

mining the star formation story of the cloud and the properties outflow nearb y another YSO driving an a lmost spherical mass 

of the nascent stars oflow, intermediate, and high mass, forming ejection (e.g., Torrelles et al. |200JJ |2003i Zapata et al.|2008), 

within it. Regarding the star formation story, there is increas- indicating that the initial conditions in the cloud may not be 

ing evidence that a single cloud can undergo different episodes the onl y a § ent determining the formation and evolution of the 

of star formation, as suggested by studies of deeply embedded mtermediate/high-mass members in a forming cluster. While 

clusters observed at spatial scales similar to the cluster member most of these cases have been found at vel T hl 8 h s P atlal scales 

separation (~ 5000 AU). In these studies, the intermediate/high- (~ 1 AU )> a broad observational base in other regions and at 

mass young stellar objects (YSOs) in the forming cluster seem othe r spatial scales is required to properly understand the role of 

to be in different evolutionary stages (judging from the peak of initial conditions of the cloud in the cluster formation process 

and ultimately in the formation of intermediate/high-mass stars. 

* The fits files for Figs.[I]|2and|4]are also available in electronic for- hi this paper we show a high angular resolution study of 

mat at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) a deeply embedded cluster whose intermediate-mass protostars 
or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-i-A/ show differences not only in their spectral energy distributions, 
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ABSTRACT 

Aims. Our main aim is to study the influence of the initial conditions of a cloud in the intermediate/high-mass star formation process. 
Methods. We observed with the VLA, PdBI, and SMA the centimeter and millimeter continuum, N2H + (1-0), and CO (2-1) emission 
associated with a dusty cloud harboring a nascent cluster with intermediate-mass protostars. 

Results. At centimeter wavelengths we found a strong source, tracing a UCH n region, at the eastern edge of the dusty cloud, with 
a shell-like structure, and with the near-infrared counterpart falling in the center of the shell. This is presumably the most massive 
source of the forming cluster. About 15" to the west of the UCH n region and well embedded in the dusty cloud, we detected a strong 
millimeter source, MM1, associated with centimeter and near-infrared emission. MM1 seems to be driving a prominent high- velocity 
CO bipolar outflow elongated in the northeast-southwest direction, and is embedded in a ridge of dense gas traced by N 2 H + , elongated 
roughly in the same direction as the outflow. We estimated that MM1 is an intermediate-mass source in the Class 0/1 phase. About 15" 
to the south of MM 1 , and still more deeply embedded in the dusty cloud, we detected a compact millimeter source, MM2, with neither 
centimeter nor near-infrared emission, but with water maser emission. MM2 is associated with a clump of N 2 H + , whose kinematics 
reveal a clear velocity gradient and additionally we found signposts of infall motions. MM2, being deeply embedded within the dusty 
cloud, with an associated water maser but no hints of CO outflow emission, is an intriguing object, presumably of intermediate mass. 
Conclusions. The UCHn region is found at the border of a dusty cloud which is currently undergoing active star formation. Two 
intermediate-mass protostars in the dusty cloud seem to have formed after the UCHn region and have different properties related to 
the outflow phenomenon. Thus, a single cloud with similar dust emission and similar dense gas column densities seems to be forming 
objects with different properties, suggesting that the initial conditions in the cloud are not determining all the star formation process. 

Key words. Stars: formation - ISM: dust - ISM: H n regions - ISM: individual objects: IRAS 001 17+6412 - Radio continuum: ISM 

1 . Introduction their spectral energy distributions: Beuther et al. 2007; Palau et 

al. l2007al b: Leurini et al. EOOTl Williams et al. 120091 . 
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Table 1. Main parameters of the VLA, PdBI, and SMA observations 
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Phase Calibrators 
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Beam 


PA. 


cont. rms 




(cm) 


Obs 


Config 


Observation 


(Bootstrapped Flux, Jy) 


Calibrators 


(" x ") 


o 


(mjy beam -1 ) 


Ref. 


6.0 


VLA 


C 


2006 Oct 27 


0102+584 (2.897 + 0.011) 


3C48 


5.2x4.0 


-4 


0.03 


1 


3.6 


VLA 


B 


1992 Jan 6 


2230+697 (0.415 ± 0.010) 


3C48 


1.0x0.8 


-66 


0.10 


3 


3.6 


VLA 


C 


2006 Oct 27 


0102+584 (3.172 + 0.029) 


3C48 


3.1 x 2.7 


+3 


0.02 


1 


3.6 


VLA 


D 


2004 Jul 10+13 


0102+584 (1.885 + 0.003) 


3C48 


14.9 x 8.5 


+76 


0.04 


2 


1.3 


VLA 


C a 


1992 May 25 


0228+673 (1.772 ± 0.056) 


3C48 


1.3 x 1.0 


+50 


0.6 


3 


1.3 


VLA 


D 


2004 Jul 10+13 


0102+584 (2.256 + 0.015) 


3C48 


5.9x4.5 


-21 


0.07 


2 


0.7 


VLA 


D 


2007 May 12 


0102+584 (3.084 + 0.121) 


3C286 


4.0x3.0 


+64 


0.17 


1 


0.32 


PdBI 


CD 


2004 Oct 17 + Dec 7 


0212+735 (1.150) 


MWC349 


4.3 x 3.4 


+61 


0.18 


1 












2145+067 










0.12 


PdBI 


CD 


2004 Oct 17 + Dec 7 


0212+735 (0.620) 


MWC349 


1.5 x 1.2 


+58 


1.0 


1 












2145+067 










0.13 


SMA 


compact 


2007 Jun 28 


0102+584 (2.584) 


Uranus 


3.2 x 2.0 


+45 


0.8 


1 



0014+612 (1.163) 

a Data from the VLA archive. The low sensitivity of the archival dataset does not allow to obtain better maps (when we combined the C and D 
configuration data) than previous VLA-D data maps (Sanchez-Monge et al. 2008 I. 
References: (1) This work, (2) Sanchez-Monge et al. 2008, (3) VLA archival data. 



but also in the ejection phenomena associated with the YSOs. 
The region, IRAS 001 17+6412, was selected from the list of 
Molinari et al. ( 119961 1 in a search for deeply embedded clus- 
ters which are luminous (> 1000 Lq) and nearby (distance 
< 3 kpc). The IRAS source has a bolometric luminosity of 
1400 Lq and is located at a distance of 1.8 kpc (Molinari et al. 
1996). The millimeter single-dish image reported by Sanchez- 
Monge et al. (2008) shows strong emission with some sub- 
structure, tracing a dusty cloud, and the centimeter emission re- 
veals an ultra-compact H n (UCH n) region, associated with the 
brightest 2MASS source of the field at the eastern border of 
the dusty cloud. The dusty cloud is associated with an embed- 
ded cluster reported by Kumar et al. (2006]!, two H2O masers 
spots (Cesaroni et al. |T9"88l Wouterloot et al. [19931 , and CO (2- 
1) bipolar outflow emission (Zhang et al. 2005; Kim & Kurtz 
2006). All this is suggestive of the dusty cloud harboring a 
UCH n region and a deeply embedded stellar cluster forming in 
its surroundings. 

We conducted high-sensitivity radio interferometric observa- 
tions in order to study the different millimeter sources embedded 
in the dusty cloud. To properly characterize the protostars, we 
also studied the distribution of dense gas, outflow and ionized 
gas emission. In the present paper we show the first results ob- 
tained from N2H + dense gas and CO outflow emission, focusing 
mainly on the intermediate/high-mass content of the protoclus- 
ter. 



2. Observations 

2. 1 . Very Large Array observations 

IRAS 00117+6412 was observed with the Very Large Array 
(VLA0) at 6 and 3.6 cm with 5-8 EVLA antennae in the ar- 
ray. The phase center of these observations was a(J2000) = 
00 h 14 m 27.725, and <5(J2000) = +64°28'46'.'171. The integration 



1 The Very Large Array (VLA) is operated by the National Radio 
Astronomy Observatory (NRAO), a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 



time was about 45 minutes at both wavelengths. Absolute flux 
calibration was achieved by observing 3C48, with an adopted 
flux density of 5.48 Jy at 6 cm, and 3.15 Jy at 3.6 cm. The 
data reduction followed the VLA standard guidelines for cali- 
bration and imaging, using the NRAO package AIPS. Images 
were performed using different weightings. The robust param- 
eter of Briggs (1995 ) was set equal to 1 and 3 (almost natural) 
respectively at 6 and 3.6 cm. The configuration of the array, ob- 
serving dates, phase calibrator and bootstrapped fluxes, synthe- 
sized beams and the continuum rms noise achieved are listed in 
Table 1. 

The observations at 7 mm (together with NH3 and SiO 
molecular line observations, Busquet et al. in prep.) were car- 
ried out with 9 EVLA antennae in the array. In order to mini- 
mize the effects of atmospheric fluctuations, we used the tech- 
nique of fast switching (Carilli & Holdaway |1997t between the 
source and the phase calibrator over a cycle of 120 seconds, with 
80 seconds spent on the target and 40 seconds on the calibrator. 
The on-source integration time was about 1.8 hours. Absolute 
flux and phase calibrations were achieved by observing 3C286 
(1.45 Jy) and 0102+584 (see Table 1). Data reduction was per- 
formed following the VLA guidelines for the calibration of high 
frequency data, using the NRAO package AIPS. The image was 
constructed using natural weighting and tapering the Mv-data at 
50 kA to increase the signal-to-noise ratio. 

Searching the archive of the VLA, we found B-array data 
at 3.6 cm and C-array data at 1.3 cm observed in 1992 (project 
AC295). The data were reduced with the standard AIPS proce- 
dures. 

2.2. Plateau de Bure Interferometer observations 

The continuum emission at 3 and 1.2 mm was observed simul- 
taneously together with the N 2 H + (1-0), and CS (5^4) molecu- 
lar transitions with the Plateau de Bure Interferometer (PdB0). 
The observations were carried out during 2004 October 17 and 



2 The Plateau de Bure Interferometer (PdBI) is operated by the 
Institut de Radioastronomie Millimetrique (IRAM), which is supported 
by INSU/CNRS (France), MPG (Germany), and IGN (Spain). 
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Fig. 1. IRAS 00117+6412 continuum maps, (a): Grey-scale: 2MASS /f s -band infrared emission. Contours: VLA-C 6 cm con- 
tinuum emission. Levels are -3, and 3 to 33 in steps of 6, times 0.03 mJy beam -1 , (b): VLA-CD 3.6 cm continuum emission. 
Levels are -3, and 3 to 30 in steps of 3, times 0.03 mJy beam -1 , (c): VLA-D 1.3 cm continuum emission. Levels are -3, 3, 
and 6 times 0.07 mJy beam -1 (Sanchez-Monge et al. 120081 1. (d): VLA-D 0.7 cm continuum emission. Levels are -3, 3, and 6 
times 0.17 mJy beam -1 , (e): Black contours: PdBI 3.2 mm continuum emission. Levels are -3, and 3 to 21 in steps of 3, times 
0.19 mJy beam -1 . Grey contours: IRAM30m continuum emission tracing the 'dusty cloud', from Sanchez-Monge et al. ([2008 ). 
Levels are 4, 6, 8, 10, and 12 times the rms of the map, 8 mJy beam -1 , (f): SMA 1.3 mm continuum emission. Levels are -3, and 3 
to 69 in steps of 6, times 0.62 mJy beam -1 . Green triangles indicate the (approximate) position of the H2O masers in the region (see 
Cesaroni et al. 1988). (g): VLA-BCD 3.6 cm continuum emission. Levels are 3 to 8 in steps of 1, times 0.05 mJy beam -1 , (h): PdBI 
1.2 mm continuum emission. Levels are -3, 3 to 1 1 in steps of 1, times 1.1 mJy beam -1 . In all panels, five-point red stars indicate 
the position of 2MASS sources associated with the main centimeter and millimeter sources. Blue dots indicate the position of the 
UCH11 region, MM 1 -main, and MM2, in order of descending declination. Synthesized beams are shown at the bottom-left corner 
of each panel, and correspond to those indicated in Table [2] 



December 7, with the array in the D (5 antennae) and C (6 
antennae) configurations, respectively. The phase center was 
o*(J2000) = 00 h 14 m 25 s .8, and 5(J2000) = +64°28'43'/00, and 
the projected baselines ranged from 24.0 to 229.0 m. The sys- 
tem temperatures were ~ 300 K for the receiver at 3.2 mm, and 
~ 2000 K for the receiver at 1 .2 mm for both days. Atmospheric 
phase correction based on the 1.2 mm total power was applied. 
The receiver at 3.2 mm was tuned to 93. 17378 GHz (lower side- 
band) to cover the N2FL (1-0) transition, for which we used a 
correlator unit of 20 MHz of bandwidth and 512 spectral chan- 
nels, providing a spectral resolution of 0.04 MHz (0. 13 km s -1 ). 



The receiver at 1.2 mm was tuned to 244.93556 GHz (upper 
sideband), covering the CS (5—4) transition, and the CH3OH 5q- 
4o A+ and 5 1 — 4- 1 E transitions. The CS and CH3OH line emis- 
sion at 1 .2 mm will be presented in a subsequent paper (Busquet 
et al., in prep.). For the continuum measurements we used two 
correlator units of 320 MHz in each band for both receivers. The 
FWHM of the primary beam was ~54" at 3.2 mm, and ~22" at 
1.2 mm. 

Bandpass calibration was performed by observing 3C454.3 
for October 17 and 2145+067 for December 7. We used the 
source 0212+735 to calibrate the phases and amplitudes of the 
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Table 2. Multiwavelength results for the intermediate-mass YSOs in the star-forming region IRAS 001 17+6412 
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0{J IVVU.U) — 


+04 Zo 40. 


















6.0 


VLA-C 


5.2 x 4.0 


-4 


30 


1.21 ± 0.03 


1.72 ± 0.08 


3.1 


± 0.5 x 2.7 ± 0.4 


95 ± 


50 


3.6 


VLA-CD 


3.8x3.7 


+45 


30 


1.03 ±0.03 


1.88 ±0.09 


3.7 


±0.4x3.1 ±0.4 


10 ± 


60 


3.6 


VLA-BCD 


1.9 x 1.7 


-70 


50 


0.43 ± 0.05 


1.62 ±0.16 


3.9 


± 0.4 x 3.7 ± 0.4 


loU ± 


An 
OU 


1.3 C 


VLA-D 


5.9x4.5 


-21 


70 


0.59 ± 0.07 


1.7 ±0.2 


5.9 


± 1.8x5 ±2 


YISj ± 


^n 
jU 


0.7 


VLA-D 


4.0x3.0 


+64 


170 


1.35 ±0.17 


1.4 ±0.3 


2.7 


± 1.2x0.0 


i cc\ _i_ 
IjU ± 


JU 


0.32 


PdBI-CD 


4.3 x 3.4 


+61 


180 




< 1 










0.13 


SMA-comp 


3.2x2.0 


+45 


600 




< 4 
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PdBI-CD 
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0(J ZUUU.U) = 


i A /I °OQ' A 1" "7b 
+ 04 ZO 4j>. / 


















6.0 


VLA-C 


5.2 x 4.0 


-4 


30 




< 0.12 










3.6 


VLA-CD 


3.8x3.7 


+45 


30 


0.14 ±0.03 


0.17 ±0.07 


6.4 


±3.0x3±2 


120 + 


30 


3.6 


VLA-BCD 


1.9 x 1.7 


-70 


50 


0.21 ±0.05 


0.15 ±0.06 


4.0 


± 1.5 x 2.2 ± 1.5 


i An -l 
10U ± 


jU 


1.3 


VLA-D 


5.9x4.5 


-21 


70 




< 0.27 










0.7 


VLA-D 


4.0x3.0 


+64 


170 




<0.69 










0.32 


PdBI-CD 


4.3 x 3.4 


+61 


180 


4.16±0.19 


5.1 ±0.9 


3.2 


± 0.4 x 0.0 


IJ ± 


1U 


0.13 


SMA-comp 


3.2x2.0 


+45 


600 


51 ± 1 


71 ± 14 


2.4 


±0.3x 1.0 ±0.2 


^n -l 
DU ± 


i n 
1U 


0.12 d 


PdBI-CD 


1.5 x 1.2 


+58 
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21 ± 1 
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±0.4x 1.6 ±0.4 
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6.0 


VLA-C 


5.2 x 4.0 


-4 


30 




< 0.12 










3.6 


VLA-CD 


3.8x3.7 


+45 


30 




< 0.12 










3.6 


VLA-BCD 


1.9 x 1.7 


-70 


50 




< 0.21 










1.3 


VLA-D 


5.9x4.5 


-21 


70 




< 0.27 










0.7 


VLA-D 


4.0x3.0 


+64 


170 




<0.69 










0.32 


PdBI-CD 


4.3 x 3.4 


+61 


180 


2.8 ± 0.2 


3.2 ±0.6 


1.6 


± 1.0 x 1.0 ± 1.5 


160 ± 


50 


0.13 


SMA-comp 


3.2x2.0 


+45 


600 


19 ± 1 


24 ±5 


2.0 


± 0.5 x 0.0 


170 ± 


20 


0.12 d 


PdBI-CD 


1.5 x 1.2 


+58 


1000 


20 ±5 


50 ± 13 


3.5 


± 1.0 x 1.1 ±0.8 


140 ± 


30 



a Primary beam corrected (mainly affecting the millimeter data). In case of non-detection at one of the frequencies, an upper limit of 4cr was 

used. Error in intensity is ltr. Error in flux density has been calculated as yj(cr V^urcc/^bcam) 2 + (ffiux-scaie) 2 , where cr is the rms noise level 

of the map, 9 somce and $beam are the size of the source and the beam respectively, and cr flux _ sca]e is the error in the flux scale, which takes into 

account the uncertainty on the calibration applied to the flux density of the source (S v x %uncertamty)- 
b Coordinates for the UCHn region, MM1, and MM2 were measured in the VLA 3.6 cm, PdBI 1.2 mm, and SMA 1.3 mm images (Fig. [TJ, 

respectively. Note that the coordinates of MM1 correspond to the coordinates of MMl-main, listed in Table[3] 
c Data from Sanchez-Monge et al. ( 120081 ). Deconvolved size and P.A. correspond to the component VLA 3 from Sanchez-Monge et al. ( 2008 1, 

which is the component associated with the easternmost 2MASS source shown in Fig.QJ. 
d 1.2 mm PdBI continuum map done with natural weighting; for flux densities and coordinates for the different subcondensations detected in 

the uniform-weighted map, see Table[3] 



antennae for both days. The rms noise of the phases was < 50° 
for the data at 3.2 mm and < 60° at 1.2 mm. The absolute 
flux density scale, calibrated with MWC349 on October 17 
and 2145+067 on December 7, has an estimated uncertainty of 
around 15%. Data were calibrated using CLIC and imaged with 
MAPPING, both part of the GILDAS software package. Imaging 
of the 3 and 1 .2 mm emission was performed using natural and 
uniform weighting, respectively. See Table 1 for details on the 
synthesized beams and the rms noises. 

2.3. Submillimeter Array observations 

The Submillimeter Array (SMAEl Ho et al. 120041 1 in the com- 
pact configuration was used to observe the 1.3 mm continuum 

3 The SMA is a joint project between the Smithsonian Astrophysical 
Observatory and the Academia Sinica Institute of Astronomy and 



emission and the 12 CO(2-l) molecular transition line (centered 
at 230.538 GHz, upper sideband) on 2007 June 28. The phase 
center of the observations was ff(J2000) = 00 h 14 m 25 s .80, and 
c5(J2000.0) = +64°28'43'.'0, and the projected baselines ranged 
from 9 to 78 kA (12-101 m). System temperatures ranged be- 
tween 80 and 200 K. The zenith opacities, measured with the 
NRAO tipping radiometer located at the Caltech Submillimeter 
Observatory, were good during the track, with r(225 GHz)~ 
0.08. The correlator, with a bandwidth of 1.968 GHz, was set 
to the standard mode, which provided a spectral resolution of 
0.8125 MHz (or 1.06 km s~' per channel) across the full band- 
width. The FWHM of the primary beam at 230 GHz was ~ 56". 
The flagging and calibration of the data were done with the MIR- 



Astrophysics, and is funded by the Smithsonian Institution and the 
Academia Sinica. 
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Table 3. Parameters of the 1.2 mm subcondensations associated with MM1 from the uniform-weighted PdBI map a 



5 





Position 


^peakb 


5/ 


Deconv. Size c 


P.A. C 


M env d 


Source 


a(J2000) 


(5(J2000) 


(mjy beam -1 ) 


(mJy) 


(" X ") 


O 


(M Q ) 


MMl-main e 


00 14 26.051 


64 28 43.69 


12.8 + 1.1 


39 + 6 


3.2+1.8x0.3 + 0.3 


-89+ 5 


1.46 












2.1 ±0.9 x 1.0 + 0.3 


+3 + 13 




MM1-S 


00 14 26.065 


64 28 42.26 


4.7 ± 1.1 


5.4+1.4 


1.6 + 0.5x0.8 + 0.2 


-30 ± 30 


0.20 


MM1-SW 


00 14 25.861 


64 28 42.25 


5.9 ± 1.1 


6.4+1.5 


1.3 + 0.4x0.2 + 0.1 


+85+ 9 


0.24 



a Synthesized beam of the uniform-weighted map is 0796 x 0777 with P.A.= +66°. The rms noise level of the map is 1.11 mJy beam -1 . 
b Primary beam corrected. Error in intensity and flux density as in Table[2] 

c Deconvolved sizes and P.A. obtained from a fit with two Gaussian for MMl-main, and one single Gaussian for MM1-S and MM1-SW. 

d Masses derived assuming a dust mass opacity coeffecient at 1.2 mm of 0.9 gcirr 1 (agglomerated grains with thin ice mantles in cores of 

densities ~ 10 6 cm 4 , Ossenkopf & Henning 1994), and a dust temperature of 30 K. The uncertainty in the masses due to the opacity law is 

estimated to be of a factor of 2. 
e Throughout all the text, we refer to MMl-main as MM1, for simplicity. 



IDL0 package. The passband response was obtained from obser- 
vations of 3C454.3. The baseline -based calibration of the am- 
plitudes and phases was performed using the sources 0102+584 
and 0014+612. Flux calibration was set by using Uranus, and 
the uncertainty in the absolute flux density scale was ~ 20%. 
Imaging and data analysis were conducted using the standard 
procedures in MIRIAD (Sault et al. [19951 ) and AIPS (see Table 1 
for details). The continuum was obtained by averaging all the 
line-free channels of the upper sideband and the lower sideband. 

3. Results 

3.1. Centimeter continuum emission 

We detected centimeter radio continuum emission at all wave- 
lengths. In Fig. QJ we show the 6 cm continuum emission 
of the region. The field is dominated by a strong and com- 
pact source associated with the brightest infrared source in the 
field, which is the counterpart of the UCH n region detected at 
3.6 cm by Sanchez-Monge et al. (2008). Additionally, the wide 
field of the VLA at 6 cm allows us to detect two unresolved 
sources to the north and east of the UCH n region. The northern 
one, with the coordinates of a(J2000.0) = 00 h 14 m 03 s .21, and 
(5(J2000.0) = +64°32'26".5, has a primary beam corrected flux 
density of 1 . 1 + 0.2 mJy. The eastern one, with the coordinates 
of a(J2000.0) = 00 h 14 m 57 s .69, and £(J2000.0) = +64°28'49".6, 
has a primary beam corrected flux density of 0.44 ± 0.09 mJy. 
Both sources are outside the region shown in Fig.[T^. 

At 3.6 cm, we improved the angular resolution of previ- 
ous observations by a factor of five. These observations reveal 
two sources in the field: an eastern source coincident with the 
UCH ii region, and a second and fainter source located ~ 15" to 
the west. Fig.QJ) shows the resulting image after combining the 
new VLA C-configuration observations with previous VLA D- 
configuration observations (Sanchez-Monge et al. 120081 1. Since 
the combined dataset has a better Mv-coverage, we recovered 
faint structure to the south of the UCH n region similar to the 
double-source structure detected at 1.3 cm by Sanchez-Monge 
et al. d20081 l. The faint source to the west of the UCHu re- 
gion is coincident with the main 1 .2 mm peak detected with the 
IRAM 30 m telescope (Sanchez-Monge et al. 2008 ). By com- 
bining the 3.6 cm VLA-C and VLA-D datasets with archival 
VLA-B data we improved the angular resolution up to ~ 1.8". 

4 The MIR cookbook by Charlie Qi can be found at 
http://cfa-www.harvard.edu/~cqi/mircook.html 



With this angular resolution, the faint source 15" to the west of 
the UCH ii region is marginally detected, and the UCH n region 
shows a shell morphology with three main peaks and with the 
2MASS source falling at the center of the peaks (see Fig. Q~g). 

In Fig. [TJ; we show, for completeness' sake, the 1.3 cm con- 
tinuum emission map from Sanchez-Monge et al. (2008jf. At 
7 mm we detected one compact source at the position of the 
UCH ii region (Fig. [TJi). 

In Table|2]we summarize the main results of the sources de- 
tected in the region at different wavelengths. The Table gives for 
each source the coordinates, peak intensity, flux density, decon- 
volved size and position angle. 

3.2. Millimeter continuum emission 

Figure QJ shows the 3.2 mm PdBI continuum emission of the 
region. There are two compact strong sources at 3.2 mm: the 
strongest source, MM1, lying ~ 15" to the west of the UCHu 
region; and the other source, MM2, being located ~ 15" to 
the south of MM1. The peak of MM1 is associated with the 
faint 3.6 cm source shown in Fig. [TJ> and with a near-infrared 
2MASS source. MM1 is partially extended to the southwest, 
also spatially coincident with a second infrared source detected 
in the 2MASS (J00142558+6428416), hereafter 2M0014256 
(cf. Fig. QJ). Note that the /f s -2MASS image shows no infrared 
emission toward MM2 (cf. FigQJ). 

The SMA map at 1 .3 mm shows two compact sources clearly 
associated with MM1 and MM2 (Fig. QJ). The 1.2 mm maps 
from the PdBI, with an angular resolution three times better than 
the SMA images, but only covering the region of MM1 within 
the primary beam of 22", show that this source splits up into at 
least three subcondensations: a compact core, MMl-main, elon- 
gated in the east- west direction and with faint extensions towards 
the north and the south; a faint 4<x source about 1 " to the south of 
MMl-main, MM1-S; and a faint source at 5cr located ~ 2" to the 
southwest of MMl-main, MM1-SW (see Fig.[TJi and TableQ. 
We note that the PdBI 1 mm and SMA 1 mm flux densities are 
completely consistent, if we take into account that the PdBI is fil- 
tering out emission at smaller scales than the SMA. Throughout 
all the paper we will refer to MMl-main as MM1. 

In Table|2]we summarize the main results of the sources de- 
tected in the region at millimeter wavelengths (see Sect. 13.11 for 
a description of the table). For both MM1 and MM2, we esti- 
mated the mass of the dust component assuming dust tempera- 
tures of 30 K and 20 K for MM1 and MM2, respectively (as a 
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Fig.2. Channel maps of the isolated N 2 H + (1-0) hyperfme F\F = 01 -* 12 line, toward the IRAS 00117+6412 region, with a 
spectral resolution of 0. 16 km s _] . The central velocity of each channel is indicated in the upper right corner, and the velocity to take 
as reference is -44.3 km s _1 . The synthesized beam, shown in the bottom left corner of each panel, is 4'.'30 x 3'.' 39, at P.A. = 60?7. 
Contours are -6, -3, and 3 to 18 in steps of 3, times the rms noise, 0.013 Jy beam 1 . Blue dots indicate the position of the UCHu 
region, MM1, and MM2 as in Fig.Q] 



first approximation, since MM1 is associated with near-infrared 
emission, while MM2 is not), and a dust mass opacity coeffi- 
cient at 1.3 mm of 0.9 g cirT 1 (agglomerated grains with thin ice 
mantles in cores of densities ~ 10 6 cm -3 , Ossenkopf & Henning 
119941 ). With these assumptions, we estimated a mass for MM1 of 
-3.0 Mq, and a mass for MM2 of ~ 1 .7 Mq (both derived from 
the (SMA) flux at 1.3 mm given in Table|2]i. In order to properly 
estimate the spectral index between 3.2 mm (PdBI) and 1.3 mm 
(SMA), we made the SMA image using the same KV-range as 
the PdBI data (9-72 kA). By measuring the flux densities in the 
maps in the common Mv-range, we obtained a spectral index of 
3.1 ± 0.2 for MM1 and 2.7 ± 0.2 for MM2. Finally, we esti- 
mated the fraction of missing flux resolved out by the SMA at 
1.3 mm by comparing the total flux of MM1 plus MM2 with the 
flux density measured with the IRAM30m Telescope of 1.2 Jy 



(Sanchez-Monge et al. 2008), and found that the fraction of flux 
filtered out by the SMA is 92%. 

3.3. N 2 H + emission 

Figure [2] shows the N 2 H + (1-0) channel maps of the hyperfme 
F\F = 01 — > 12 line (hereafter, the 'isolated' hyperfme) toward 
IRAS 001 17+6412. Close to the systemic velocit>0 and span- 
ning ~ 2 km s~', the emission appears as a ridge at the MM1 



5 The syste mic ve locity for IRAS001 17+6412 is -36.3 km s~' 
(Molinari et al. [19961 Zhang et al. l2005l . and the N 2 H + (1-0) hyperfme 
adopted as reference line in this work was F\F = 23 — > 12. However, 
since the hyperfme adopted as reference is blended with two other hy- 
perfines, we will make our analysis with the hyperfme which is isolated, 
the F\F = 01 — > 12 line, which is at 8 km s~' with respect to the line 
adopted as reference. Therefore, the reference velocity (corresponding 
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Fig. 3. a) Zero-order moment (integrated intensity) for the 
hyperfine F r F = 01 -> 12 line of N 2 H + (l-0) toward 
IRAS 00117+6412. Contour levels range from 1 to 91% of the 
peak intensity, 0.357 Jy beam -1 km s _1 , increasing in steps of 
15%. Grey-scale: 2MASS /T s -band infrared emission, b) First- 
order moment (velocity) for the hyperfine F\F = 01 — > 12 line 
of N2H + (l-0) (wedge units are km s _1 ). Note that the velocity 
to take as reference is -44.3 km s _1 . Note also that the cross 
corresponds to the NE position of the MM1 ridge (4.0", 2.5" 
from phase center), where a spectrum is shown in Fig. [6] and 
is used to perform the pv-plot at PA=130° shown in Fig. [7^. 
c) Second-order moment (velocity dispersion) for the hyperfine 
F\F = 01 -> 12 of N 2 H + (l-0) (wedge units are km s" 1 ). Note 
that the velocity dispersion must be multiplied by the factor 4.1. kinematics: moments and pv-plots 



MM2 clump), the most blueshifted velocities show some faint 
emission ~20" towards the west of MM2 (western clump). The 
zero-order moment map integrated for the isolated hyperfine of 
N2H + (l-0) is presented in Fig. [3^, overlaid on a /T s -band image 
of 2MASS. The MM1 ridge has a length of 20" and the crest is 
coinciding well with the peak and elongation of the 3.2 mm con- 
tinuum emission (Fig. QJ). As for the MM2 clump, it has a size of 
~ 10", with the peak position also coinciding with the millime- 
ter peak, and is elongated in the southeast-northwest direction 
(PA~ 132°). Finally, it is worth noting that the MM2 clump and 
the western clump fall at a region which is completely dark in 
the near-infrared, as seen in the /f s -2MASS image of Fig.[3k. 

3.4. CO (2-1) emission 

Channel maps of the CO (2-1) emission are displayed in Fig. [4] 
The CO (2-1) emission appears spanning a wide range of ve- 
locities, from -70 up to -14 km s (systemic velocity at 
-36.3 km s~'). Blueshifted emission appears in different clumps 
toward the north and northeast of MM1, while redshifted emis- 
sion appears toward the southwest of MM1, suggesting a bipolar 
structure centered on MM1. 

The spectrum of the emission integrated over all the region is 
shown in Fig. [5] together with a preliminary spectrum of I3 CO 
(Busquet et al., in prep.). The CO spectrum shows a dip from 
-6 up to +6 km s _1 with respect to the systemic velocity, which 
could be due to self-absorption by cold foreground gas or opac- 
ity effects, as the dip is coincident with the peak of the I3 CO 
line. We note that the dip could be partially produced by the 
missing short-spacing information in the interferometer data as 
well. In order to make a rough estimate of the fraction of flux 
filtered out by the SMA (due to the lack of uv sampling at spac- 
ings smaller than 9 kA), we compared the single-dish spectrum 
(from Zhang et al. 2005) at certain velocities with the SMA spec- 
trum (extracted from the channel maps of Fig. |4]convolved to the 
single-dish beam, of 30") at the position of the IRAS source. The 
missing flux is around 99% for the systemic velocities, indicat- 
ing that at these velocities the CO emission mainly comes from 
structures much larger than ~ 10", which is the largest angu- 
lar scale observable by an interferometer whose shortest base- 
line is ~ 9 k/1 (see Appendix for details). On the other hand, 
the missing flux decreases as one moves to higher velocities. 
For example, the fraction of missing flux goes down to 50% at 
-7 and 5 km s _1 (with respect to the systemic velocity) for the 
blueshifted and redshifted sides, respectively; and the SMA re- 
covers all the single-dish flux at -13 and 7 km s _1 with respect 
to the systemic velocity. This indicates that the high-velocity 
CO (2-1) emission has a characteristic source size < 10". 



4. Analysis 



2 V2Zn2 =i 2.35 to have linewidths. In all panels, black/white 
dots indicate the position of the UCH 11 region, MM1, and MM2 
as in Fig.Q] 



position (labeled as the MM1 ridge) and is elongated roughly 
in the east-west direction (PA~ 58°). While at redshifted veloc- 
ities the emission is only associated with MM2 (hereafter, the 



to the systemic velocity) in the analysis of the N2H + emission from the 
isolated hyperfine is -36.3 - 8.0 = -44.3 km s -1 . 



In order to study the kinematics of the dense gas as traced by 
N2H + , we computed the first-order moment for the isolated hy- 
perfine of N2H + (l-0), which does not suffer from blending with 
other hyperfines (Fig. [3})). The first-order moment map shows 
a velocity gradient in the MM1 ridge, spanning velocities of 
-43.5 to -45.0 km s , centered approximately on MM1, and 
perpendicular to the elongation of the ridge. Another clear ve- 
locity gradient can be seen toward MM2 at red velocities, from 
-42 to -44 km s _1 . Note that both velocity gradients observed 
toward MM1 and MM2 have the blueshifted velocities toward 
the northwest and the redshifted velocities toward the southeast. 
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Fig. 4. CO(2-l) channel maps of the IRAS 00117+6412 region, averaged over 3.1 km s _1 wide velocity intervals. The central 
velocity of each channel is indicated in the upper left corner, and the systemic velocity is -36.3 km s _1 . Symbols are the same as in 
Fig-Ill The synthesized beam, shown in the bottom left corner of each panel, is 2'.'92 x 1'.'78, at P.A. = 44°.7. Contours are -4, and 
4 to 40 in steps of 6, times the rms noise, 0.1 Jy beam -1 . 



In Fig. [3]: we show the second-order moment map of the region, 
tracing the velocity dispersion, which is largest toward MM1 and 
MM2, and in general along the MM1 ridge, reaching (linewidth) 
values of 0.9-1 .4 km s , much larger than the thermal linewidth 
for N2H + , of around 0.22 km s at 30 K. The western clump is 
the more quiescent emission in the region, with linewidths of 
up to 0.6 km s _1 . In addition, the first-order moment reveals an 
abrupt change in velocity in the western clump (Fig. 03), which 
could be indicative of the clump being kinematically decoupled 
from the MM1 ridge and the MM2 clump. 



In Fig.|6]we show the N2H + spectra for MM1, for a position 
in the northeastern (NE) side of the MM1 ridge (see positions 
in Fig. [3p), and for MM2. Note that the lines are broad, as can 
be seen in the isolated hyperfine line (at ~ -44 km s _1 ), and 



suggestive of a double velocity component. See Sects. I4T21 and 
|4.3| for a further analysis and discussion on these spectra. 



We additionally built position-velocity (pv) cuts to further 
study the kinematics of the region. In Fig. [7] we show two cuts 
in the direction of the velocity gradient seen in the MM1 ridge 
(at PA~ 130°, for MM1 and a position to the NE of MM1, see 
positions in Fig. and one cut along the elongation of the 
ridge (PA=40°). The three cuts show emission spanning from 
~ -43 to ~ -45 km s _1 , with one compact component at the 
blueshifted velocities of (-44.2)-(-45.0) km s , and a broad 
red component at ~ -43.5 km s _1 . The emission from the NE 
position and from MM1 follow the same general trend of a blue 
component with a small velocity gradient and a red broad com- 
ponent with no gradient (better seen in the cut of the NE posi- 
tion, Fig. [7^). For the perpendicular cut (PA=40°, centered on 
MM1, Fig. Efc), we also find a velocity gradient with the red- 
shifted velocities at negative values, or to the southwest of the 
cut, and the blueshifted toward the northeast. Regarding MM2, 
the velocity gradient is also measured at a PA=130°, coincident 
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Velocity c 


Linewidth c 




T 

1 ex 


W(N,H + ) 


Object 


Line a 


(Kkms" 1 ) 


(km s" 1 ) 


(kms- 1 ) 


T b 


(K) 


(cm" 2 ) 


MM1 


1 


0.062 ± 0.003 


-36.7 


0.66 


0.20 ± 0.02 


3.1 


6.7 x 10" 




2 


0.060 ±0.001 


-35.9 


1.40 


<0.1 


3.4 


8.5 x 10" 


MM1-NE 


1 


0.072 ± 0.003 


-36.7 ±0.1 


0.68 ± 0.02 


0.15 ±0.02 


3.3 


5.8 x 10" 




2 


0.038 ± 0.001 


-35.6 ±0.1 


0.98 ± 0.04 


< 0.1 


3.2 


5.3 x 10" 


MM2 




0.031 + 0.001 


-34.9 ±0.1 


1.49 ±0.05 


< 0.1 


3.1 


7.6 x 10" 



11 For MM1 and the NE position of the MM1 ridge, the hyperfine structure is fitted with two velocity components. For MM2 the fit is made with 
one single velocity component. 

b A = f[J v (T cx ) - J v (Tb g )]. r m is the optical depth of the main hyperfine, adopted to be the F\ F = 23 — > 12 line. The otpical depth of the isolated 

line can be obtained by multiplying r m by the factor 3/7. 
c When errors are not given, the velocity and linewidth have been fixed from a fit of two Gaussian to the isolated line. When errors are given, 

the fit of two Gaussian to the isolated line was used to provide initial values for the velocity and linewidth of the hyperfine fit. 
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-30.0 km s _1 up to -8.0 km s _1 , and blue wing: -72.0 km s _I up 
to -42.5 km s"" 1 ) used to estimate the parameters of the outflow. 



with the elongation of the MM2 clump (Fig. [3J)), and is well 
seen in the pv-plot at PA=130°, with a double peak (Fig. [8} top). 
However, the perpendicular cut does not show any gradient, al- 
though the double peak is also there. We refer the reader to the 
subsequent sections for the interpretation of the kinematics of 
the N 2 FT dense gas toward MM1 and MM2. 

4.2. N 2 H + kinematics: MM1 ridge 

In order to fit the hyperfine spectra toward MM1 and the NE 
position of the ridge (spectra shown in Fig. |6]), we assumed 
two velocity components. We initially fitted two Gaussian to the 



isolated line to estimate the initial values of the velocities and 
linewidths for the two velocity components. With these initial 
values, we fitted the hyperfine spectra with two velocity com- 
ponents for the NE position and MM1. The results are listed in 
Tableland shown in Fig. [6] The fits can reproduce the data rea- 
sonably well. 

For MM1, the derived two velocity components at 
-36.7 km s , with a linewidth of ~ 0.7 km s _1 , and 
-35.9 km s _1 , with a linewidth of ~ 1.4 km s _1 , were already 
seen in the pv-plots of MM1 at PA=130° shown in Fig.|7](note 
that in Fig. |7]we show the isolated line, which was not taken as 
the reference line, and the component fitted at -36.7 km s _1 cor- 
responds to -44.7 km s _1 , while the component at -35.9 km s _1 
corresponds to -43.9 km s _1 in the pv-plot). For the NE posi- 
tion of the ridge, the two velocity components are found again 
around -36.7 and -35.6 km s _1 , and the redshifted component 
has a broader linewidth than the blueshifted component, as in 
the MM1 spectrum. We note that the red component is optically 
thinner than the main component. This red component, which 
shows no velocity gradient, could be tracing a filament inter- 
cepting the line of sight maybe connecting the MM1 ridge and 
MM2, as the velocity of this red component is similar to the ve- 
locity of MM2. Finally, we showed in Sect. @7T](Fig. EDO that the 
blue component (corresponding to the strongest component in 
the spectra) has a small velocity gradient which could be tracing 
the global rotation of the MM1 ridge. 

In Table |4] we also list the excitation temperature derived 
from the fits and the N2H + column density calculated assum- 
ing that the filling factor is ~ 1 and following Caselli et al. 
(2002b). We note that the excitation temperature of ~ 3 K, and 
the N2H + column densities are similar to the values found in 
low-mass star-forming regions (e. g., Caselli et al. |2002al Chen 
et al. 120071 Kirk et al. 120091 , and slightly smaller than the values 
found in massive star-forming regions observed with interferom- 
eters (e.g., Palau et al. |2007al Beuther & Henning |2009l . The 
small excitation temperature obtained is suggestive of either low 
density and/or cold gas. Note however that our values could be 
affected by the missing flux problem caused by the lack of short 
Mv-spacings in the interferometric data (for the PdB, the largest 
angular scale detectable is ~ 11", see Appendix, while in Palau 
et al. (|2007al > and Beuther & Henning (2009), the largest angu- 
lar scale was 20-30", which means that our PdB data are more 
affected by the missing flux problem). 
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Table 5. Output parameters of the position-velocity plot modeled emission for MM2 



Parameter 


Value 


Fixed/free parameter and comments 


Disk central velocity (km s ^) 


-43.0 


Fixed from observed PA— 40° pv-plot 


Disk outer radius (", adopted as reference radius) 


+5 


Fixed from emission size 


Disk inner radius (") 


+ 1 


Fixed: upper limit from SMA data (see main text) 


Intensity power law index 


-1 


Fixed following other studies of disk-like structures (see main text) 


Infall velocity power law index 


-0.5 


Fixed, free-fall 


Rotation velocity power law index 


+ 1 


Fixed, solid rigid rotation 


Reference infall velocity (km s -1 ) 
Reference rotation velocity (km s ) 


0.17 ±0.02 
0.6 ±0.2 


Free: fitted in the PA=40° pv-plot 

Free: fitted in the PA=130° pv-plot, after fitting the ref. infall vel. 



Table 6. Physical parameters of the outflow driven by MM1 





*dyn 


size 




M out a 


M" 






r 


77 a 


1 a 

*-tnech 


Lobe 


(yr) 


(arcsec) 


(cm- 2 ) 


(Mq) 


(M Q yr~ 


') 


(Mq km s _1 ) 


(Mq km s" 1 yr -1 ) 


(erg) 




Red 


3600 


12 x 9 


3.1 x 10 16 


0.05 


1.5 x 10" 


-5 


1.2 


3.3 x 10~ 4 


2.6 x 10 44 


0.46 


Blue 


3300 


14 x 7 


2.9 x 10 16 


0.05 


1.4 x 10" 


-5 


1.4 


4.2 x 10~ 4 


4.2 x 10 44 


0.85 


All 


3500 


26x9 


6.0 x 10 16 


0.10 


2.9 x 10" 


-5 


2.6 


7.5 x 10~ 4 


6.7 x 10 44 


1.31 



a Parameters are calculated for an inclination with respect to the plane of the sky equal to 0°, and are corrected for opacity effects. A' refers to 
CO column density. 



4.3. N 2 H + kinematics: MM2 

We fitted the N 2 H + (1-0) spectrum toward MM2 with a single 
velocity component (Fig. |6}bottom). As can be seen in the fig- 
ure, there is some excess of emission at blueshifted velocities, 
which could be contamination from the MM1 ridge (which has 
an extension toward the south at blueshifted velocitites with re- 
spect to the MM2 clump, at around (-44.6)-(-44.7) km s" 1 , as 
seen in Fig. [2]). Alghough the excess of emission is quite clear, a 
fit with two velocity components for the MM2 N2H + hyperfine 
spectrum could not be well determinec0. From the parameters 
obtained from the fit, we derived the excitation temperature and 
N2H + column density as in the previous section, and the values 
found for MM2 are very similar to the values found for MM1. 

In order to further constrain the N2H + kinematics toward 
MM2, we approached the flattened structure seen in N2H + to 
a model of a spatially infinitely thin disk seen edge-on, with a 
power-law intensity as function of radius, consisting of a su- 
perposition of optically thin rings undergoing infall and rotation 
(also described as power-laws). We computed the synthetic pv- 
plots along the projected major and minor axes of the disk, with 
angular and spectral resolutions of 3.4" and 0.35 km s , respec- 
tively. The modeled emission is shown in Fig.[8}bottom, and the 
adopted/fitted parameters are listed in Table [5] 

The disk structure is modeled with an inner and outer radius. 
The outer radius is constrained from the N2H + emission in the 
observed pv-plot, and is also taken as the reference radius for 
the rotation and infall power-laws. The intensity power-law in- 
dex was adopted to be -1, as used in other models of disk-like 
structures around intermediate/high-mass YSOs (e. g., Beltran 
et al. 120041 ) . For the infall we assumed a free-fall infall veloc- 



6 This is because the two-velocity-components hyperfine fit is 
strongly dependent on the adopted initial values. For the case of MM1 
and the NE position of the MM1 ridge, this problem was solved by fit- 
ting two Gaussian to the line which is isolated, for which both velocity 
components are well detected. However, in the case of MM2, the ve- 
locity component at -44.5 km s -1 is barely detected in the isolated line 
(see Fig.[6fbottom), hindering the estimate of the initial values. 



ity, and for the rotation law we assumed a solid rigid rotation, 
as recent studies show that these initial conditions are able to 
account for most observational properties of star-forming cores 
(e. g Tscharnuter et al. [20091 Walch et al. 2009; Zhilkin et al. 
2009). The free parameters are the reference infall velocity and 
the reference rotation velocity. 

We constrained the upper limit of the inner radius from the 
SMA map at 1 mm, which has an angular resolution of ~ 2.5", 
and hence any inner gap must have a radius < 1 .2". On the other 
hand, the N2H + data further constrain the inner radius to 0.8", as 
a smaller inner radius clearly cannot reproduce the double peak 
in the PA= 1 30° pv-plot. We conclude that the inner radius covers 
the range of 0.8-1.2", and adopted the mean value of ~ 1". We 
derived the reference infall velocity from the PA=40° pv-plot 
(only affected by infall), of 0.17 ± 0.02 km s" 1 . We note that 
to include infall in the model is necessary in order to reproduce 
the double peak seen in the pv-plot at PA =40°. Finally, we fitted 
the rotation velocity from the PA=130° pv-plot (affected by both 
infall and rotation). We note that in the PA=130° pv-plot the 

velocity component at 43.5 could be contaminated by the 

MM1 ridge, and hence we focused our fit on the two central 
peaks at around -43.0 km s" 1 . The resulting reference rotation 
velocity is 0.6 + 0.2 km s _1 . We also note that if we fit the data 
adopting an inner radius of 0.8" or 1.2", the new parameters are 
similar, within the uncertainties, to the parameters derived for an 
inner radius of 1". 

From the fitted infall velocity, one can estimate the mass 
of the central object by using M = v^ nf R/(2Gsin 2 i), with i be- 
ing the inclination angle to the plane of the sky. Thus, taking 
v; n f ~ 0.2 km s _1 at a reference radius of ~ 5", and assuming 
an inclination of 90°, we obtain a mass ~> 0.2 Mq, where the 
lower limit accounts for the inclination assumption. We note that 
the true mass however is likely not larger than a few tenths of so- 
lar mass, because we estimated the inclination of the MM2 N2H" 1 " 
core from its deconvolved major and minor axis (8" x 1"), and 
obtained an inclination, assuming it is intrinsically circular, of 
~ 80°, which is very close to the assumed inclination. It is in- 
teresting to compare the mass derived from infall, which is in- 



Palau et al.: Three intermediate-mass YSOs in the making 



11 




N 2 H (1-0) 

pv— plots 
toward the 
MM1 ridge 



Velocity (km/s) 



Fig. 6. Top: N2FF (1-0) hyperfine spectra toward MM1 . Middle: 
idem toward the NE position of the MM1 ridge (4.0,2.5 with 
respect to the phase center; see Fig. [3p). Bottom: idem toward 
MM2. For MM1, and the NE position, we fitted the hyperfine 
structure with two velocity components, and for MM2 we used 
one single velocity component. In the middle panel, the blue ar- 
rows indicate the position of the seven hyperfines (with statis- 
tical weights from Womack et al. 119921 and frequencies from 
Caselli et al. 1995 ) for the velocity component at -36.7 km s _I 
(Table 2]). The different hyperfine lines correspond, from left 
(blue velocities) to right (red velocities), to the quantum num- 
bers: F\F = 01 -> 12 (isolated line), 21 -> 11, 23 -> 12 (refer- 
ence line), 22 11, 11 10, 12 -» 12, and 10 -> 11. 



dicative of the mass already accreted onto the protostellar core, 
with the mass derived from the millimeter continuum emission 
for MM2 (around 1.7 Mq, Sect. 13.2b . which is tracing the mass 
of the envelope and disk surrounding the central core. For Class 
sources it is typically assumed that these two masses should be 
of the same order. Thus, the obtained accreted mass smaller than 
the disk+envelope mass is suggesting that the central protostel- 
lar object in MM2 has just recently formed. 

4.4. CO kinematics: moments and pv-plots 

We constructed the zero-order moment maps of the CO emission 
integrated over three different velocity ranges: (1) systemic ve- 
locities from -6 to +6 km s _1 with respect to the systemic veloc- 
ity (Fig.|9}top), (2) moderate velocities from -6 to -15 km s _1 , 
and from +6 to +15 km s _1 with respect to the systemic velocity 
(Fig.|9}middle), and (3) high velocities from -15 km s _I to most 
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Fig. 7. Position-velocity (pv) plots of the N2I-F isolated line be- 
tween -46 and -42 km s _1 , toward positions of the MM1 ridge, 
a) Pv-plot in a cut with PA=130°, and centered on the NE posi- 
tion of the MM1 ridge (4'.'0, 2'.'5 with respect to the phase center, 
see Fig.[3p). b) Pv-plot in a cut with PA=130°, and centered on 
MM1. c) Pv-plot in a cut with PA=40°, centered on MM1. Note 
that the velocity to take as reference is -44.3 km s _1 . For all 
panels, contours are -3, -2, 2 to 10 in steps of 1, times the rms 
noise level, 0.02 Jy beam -1 . 



negative velocities, and from +15 km s _1 to most positive veloci- 
ties with respect to the systemic velocity (Fig.|9}bottom). At sys- 
temic velocities the emission appears concentrated in two main 
structures, a clump associated with the N2H + MM1 ridge, and 
a multipeaked-arch located to the north of MM 1, just bordering 
the N2H + emission. At moderate and high velocities, which are 
less affected by the missing short spacings in the interferometric 
data, the emission has a bipolar structure centered near the posi- 
tion of MM1 and is elongated in the northeast-southwest direc- 
tion (see Fig. [9} middle and bottom). At blue moderate velocities, 
the emission splits up into two main peaks: Ba, close to the posi- 
tion of MM1 and elongated in the northeast-southwest direction, 
and Bb, located ~ 10" northwards of MM1 and with a round 
shape. A third fainter clump, Be, is located ~ 10" to the north- 
east of MM1 . Regarding the redshifted emission, it appears as an 
elongated cone-like structure, with one main clump, Ra, which is 
well aligned with the Ba clump, with a position angle of ~ 45°, 
and centered ~ 7" toward the southwest of MM1. At higher ve- 
locities, the emission is found toward Ba, Bb, Ra, and toward 
a second redshifted clump, Rb, which was not very prominent 
at moderate velocities. Clumps Bb and Rb are aligned at a po- 
sition angle of ~ 20°, and centered near MM1. Note that the 
position angles of the high-velocity CO emission are almost per- 
pendicular to the direction of the N2H + velocity gradient, found 
at ~ 130°. 

We performed different pv-plots close to MM1 at PA=45°, 
and PA=20°(see Fig. [TOb . in the two directions indicated in 
Fig. [9] Both pv-plots have been convolved with a Gaussian 
(5" x 2" at PA=135° and 110°, respectively). In both plots, 
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Table 7. 2MASS sources associated with the centimeter and/or millimeter emission toward the star-forming region 
IRAS 00117+6412 





2MASS 


Position 












J-H 


H-K s 




Source 


Identification 


a(J2000) 


<S(J2000) 


J 




H 




K s 


color" 


color" 


ip b 


UCHu 


00142828 +6428471 


00 14 28.29 


64 28 47.19 


14.94 ± 


0.08 


13.28 + 


0.07 


12.14 ±0.05 


1.66 


1.14 


+0.18 


MM1 


00142616 + 6428444 


00 14 26.16 


64 28 44.44 


15.49 ± 


0.07 


13.58 ± 


0.06 


12.22 ± 0.04 


1.91 


1.37 


+0.27 


2M0014256 


00142558 + 6428416 


00 14 25.59 


64 28 41.70 


16.71 ± 


0.14 


14.35 + 


0.05 


13.27 ± 0.04 


2.36 


1.09 


-0.27 



a Uncertainties in the colors are typically ~ 0.09 mag. 

b The infrared excess is measured as the difference between the (H-K s ) color and the (H-K s ) color corresponding to a reddened main-sequence 
star (following the reddening law of Rieke & Lebofsky [T985l >, i. e., IR excess = (H - K s ) - 0.576 x (J - K s ). 
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Fig. 8. Position-velocity (pv) plot of the N2H + isolated line to- 
ward MM2 in a cut of PA=130° (top left) and PA =40° (top 
right) and the corresponding modeled emission (bottom) of a 
disk structure with infall and rotation, adopting the parameters 
shown in Table [5] for the case of an inner radius of ~ 1 " (we 
shifted the disk with respect to the MM2 position by -1" in 
the PA=40 axis). For the top panels, (grey) black contours are 
(-16%, 16%,) 32%, 48%, 64%, 80%, and 96% of the peak emis- 
sion, 0.146 Jy beam -1 for the top-left panel, and 0.148 Jy beam" 1 
for the top-right panel. For the bottom panels, contours are 16%, 
32%, 48%, 64%, 80%, and 96% of the peak emission, which is 
in arbitrary units. 



positive and negative positions correspond respectively to north- 
eastern and southwestern positions with respect to MM1. In the 
first cut (45°) we can distinguish a strong velocity gradient from 
clump Ra to clump Ba, with velocities increasing with distance, 
following a Hubble-law pattern (up to 34 km s" 1 with respect to 
the systemic velocity in the case of clump Ba). We note that in 
the cut at PA=45°, the clump Be shows moderate velocity emis- 
sion of up to 14 km s" 1 with respect to the systemic velocity. 
Regarding the pv plot at PA=20°, the redshifted emission and 
the blueshifted clump Ba also follow a Hubble-law, as in the cut 
at 45°. 

We calculated the energetics of the outflow for each blue and 
red lobe separately (for the blue lobe we included both clumps 
Ba and Bb), assuming that all the emission comes from a single 
outflow, and listed the values in Table [6] The expressions used 
to calculate the outflow CO column density, N(CO), from the 



transition J — > 7-1, and the outflow mass, M out , are given in 
Palau et al. (2007b). We adopted an opacity in the line wings of 
~ 2 (from preliminary I3 CO data), and an excitation temperature 
of ~ 7 K, estimated from the spectrum in Fig. [5] (assuming that 
CO is optically thick and adopting a line temperature of 2.7 K). 
For the red lobe we integrated from -30 to -8 km s" 1 , and for the 
blue lobe from -72 to -42.5 km s" 1 . The dynamical timescale of 
the outflow, fdyn, was derived by dividing the size of each lobe by 
the maximum velocity reached in the outflow with respect to the 
systemic velocity (28.3 km s" 1 for the red lobe, and 35.7 km s" 1 
for the blue lobe). 

4.5. IRAS, MSX and 2MASS emission 

With the aim of finding the possible infrared counterparts of the 
sources studied in this work, we searched the IRAS, MSX and 
2MASS surveys (the region has not been observed by Spitzer). 
While the IRAS position error ellipse makes it difficult to disen- 
tangle the contribution from the UCH n region and MM1, we su- 
perposed the MSX emission on our millimeter maps and found 
that the MSX emission is compact and clearly peaking at the 
UCH ii region. Thus, this is suggestive of the major part of the 
IRAS and MSX fluxes, at least up to ~ 20 //m, coming from 
the UCH ii region. However, at 60 and 100 fim the contribution 
of the UCH ii region and MM1 is not clear, and for this reason 
we refrained from building the spectral energy distributions. As 
for the 2MASS Point Source Catalog (Skrutskie et al. 120061 1. in 
Table|7]we show the 2MASS photometry for the 2MASS coun- 
terparts of the UCH ii region, MM1, and for 2M0014256. From 
the 2MASS magnitudes, we estimated the ( J - H) and (H - K s ) 
colors and measured the infrared excess as the difference be- 
tween the (H - Ks) color and the (H - Ks) color correspond- 
ing to a reddened main-sequence star (following the reddening 
law of Rieke & Lebofsky |1985l l. The 2MASS sources associated 
with the UCH ii region and MM1 show a moderate infrared ex- 
cess typical of Class II sources, while 2M0014256 seems to be 
a reddened main-sequence star, or a Class III source. A detailed 
analysis of the low-mass content of the forming cluster in IRAS 
00117+6412 will be presented in a subsequent paper (Busquet 
et al., in prep.). 

5. Discussion 

5.1. A shell-like UCH ii region 

The centimeter range (from 6 cm up to 7 mm) of the spectral 
energy distribution of the UCH n region can be fitted assuming 
free-free optically thin emission from ionized gas with a spectral 
index of -0.03 ± 0.08 (see Fig. [TT] the spectral index was cal- 
culated from the fluxes measured in the images at 6, 3.6, 1.3 cm 
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Fig. 9. Top: Contours: CO (2-1) moment-zero map integrated 
for systemic velocities (between -6 and +6 km s _1 with re- 
spect to the systemic velocity, -36.3 km s _1 ). Levels start 
at 24%, increasing in steps of 15% of the peak intensity, 
22.9 Jy beam -1 km s . Middle: Idem for moderate velocities. 
Blue (solid) contours: emission in the velocity range from -6 to 
-15 km s with respect to the systemic velocity. Red (dashed) 
contours: emission in the velocity range from +6 to +15 km s _I 
with respect to the systemic velocity. Levels start at 9%, increas- 
ing in steps of 15% of the peak intensity, 32.7 Jy beam -1 km s _I 
and 25.8 Jy beam -1 km s _I for the blue and red-shifted ve- 
locities, respectively. Bottom: Idem for high velocities. Blue 
(solid) contours: emission in the velocity range from -15 km s _I 
to most negative values with respect to the systemic velocity. 
Red (dashed) contours: emission in the velocity range from 
+ 15 km s _1 to most positive values with respect to the systemic 
velocity. Levels start at 9%, increasing in steps of 15% of the 
peak intensity, 13.0 JybeairT 1 kms -1 and 11.1 JybeairT 1 kms -1 
for the blue and red-shifted velocities, respectively. In all pan- 
els, grey scale: N2H + zero-order moment map; grey solid lines 
indicate the possible orientation (PA=45° and 20°) of the out- 
flow(s). Black/white dots indicate the position of the UCHii re- 
gion, MM1 (main, south, and southwest), and MM2. The synthe- 
sized beam, 2.' 92 x 1'.'78 at P.A. = 44!7, is shown in the bottom 
left corner of each panel. 
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Fig. 10. Top: CO (2-1) pv-plot in the northeast-southwest di- 
rection with PA=45°, at the offset position (3'.'36; 0796), with 
respect to the phase center. Bottom: CO (2-1) pv-plot in the 
north-south direction with PA=20°, at the offset position (2'.'07; 
-0704), with respect to the phase center. For both panels, con- 
tours are -12%, -6%, 6% to 96% in steps of 10% of the 
peak intensity, 3.63 Jy beam 1 for the plot at PA=45°, and 
3.59 Jy beam 1 for the plot at PA=20°. Both pv-plots have been 
convolved with a gaussian of 5'.'0 x 270, at P.A. = 135°0 and 
P.A. = 110°0 respectively. 



and 7 mm made with the wv-range shared at all the wavelengths, 
4.35-55 kA). We calculated the physical parameters of the Hn 
region at 3.6 cm (in the CD configuration) assuming the emis- 
sion is optically thin. The results are consistent with a UCH n 
region with a size of 0.03 + 0.01 pc, a brightness temperature of 
12.4 K, an emission measure of 1.0 x 10 5 cirT 6 pc, an electron 
density of 2000 cirT 3 , a mass of inonized gas of 4.2 x 10~ 4 Mq 
(estimated from the beam averaged electron density and the ob- 
served size of the source), and a flux of ionizing photons of 
3.8 x 10 44 s . These parameteres are consistent with a UCHn 
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region driven by an early-type B2 star (from Panagia [19731 1, as 
already stated by Sanchez-Monge et al. ( 2008b. 

The high angular resolution centimeter observations pre- 
sented in this work reveal that the UCHn region reported in 
Sanchez-Monge et al. ([2008) already has a shell-like structure. 
The average radius of the shell (from an average of the three 
subcondensations) is ~ l'/2, which corresponds to 2200 AU or 
0.01 pc. From this radius, one can make a rough estimate of the 
dynamical timescale, either assuming a classical expansion of 
an HII region at ~ 10 km s _1 , or the expansion of a wind-driven 
bubble (as typically assumed to explain the shell-like morpholo- 
gies of UCHn regions, e. g., Garay & Lizano [T999l l. In the first 
case (classical HII region), the dynamical timescale is ~ 1000 yr, 
while in the second case (stellar wind blown bubble) the dynam- 
ical timescale is ~ 12000 yr (assuming the stellar wind domi- 
nates over the classical expansion of the UCH n region, an ini- 
tial ambient density of ~ 10 7 cm -3 , and following Castor et al. 
119751 and Garay & Lizano 1999 yielding an expanding veloc- 
ity < 1 km s _1 ). The lifetimes derived from both assumptions 
are much smaller than ~ 1 Myr, the lifetime estimate for UCH n 
regions of ~ 1000 Lq, from the RMS survey by Mottram et al. 
2009 (in prep.). This could be indicative of the UCH n region 
in IRAS 00117+6412 having undergone in the past a period of 
strong quenching (due to a high-mass accretion rate, Walmsley 
et al. 119951 1, before expanding as a classical HII region or as a 
wind-blown bubble. Alternatively, if the classical HII region ex- 
pansion or the wind-blown bubble assumptions are correct, we 
could be witnessing the very first stages of the expansion of an 
ionized shell around a B2-type (~ 1000 Lq) YSO. However, the 
fact that we do not detect significant dust or dense gas emission 
associated with the UCH n region suggests that the object is in 
the process of disrupting its natal cloud, and most likely has a 
lifetime similar to the < 1 Myr given by Mottram et al. (2009, 
in prep.). In addition, diffuse optical emission can be seen to the 
northeast of the UCH n region in the blue and red plates of the 
Palomar Observatory Sky Survey II, supporting the fact that the 
UCH n region is not deeply embedded and is emerging from its 
natal cloud. 

From the bolometric luminosity and the spectral type of 
the UCH ii region one can estimate the mass and the evolu- 
tion time from the birthline by placing the object in a HR dia- 
gram. Following the models of Palla & Stahler (|1990l|19931>, and 
adopting an effective temperature of ~ 20500 K (Panagia l 19731 1, 
the mass of the star ionizing the UCH n region is about 5-6 Mq, 
and the star is already placed at the ZAMS, which most likely 
is reached after ~ 1 Myr for a star of this mass (e. g., Palla & 
Stahler [T9931 Bernasconi & Maeder [l996l l. This estimate agrees 
with the previous estimate of the lifetime of a UCHn region 
of ~ 1000 Lq (derived from Mottram et al. 2009, in prep.) of 
~ 1 Myr. 

5.2. MM1: multiple sources and a spectacular outflow 

In previous sections we showed that MM1 is associated with a 
2MASS source with infrared excess. Dense gas emission traced 
by N2H + and faint centimeter emission are also associated with 
MM1. These properties suggest that MM1 is most likely in 
a phase of moderate/intense accretion, possibly in the Class 
0/1 stage, with typical ages of (0.5-l)xl0 5 yr (e.g., Evans et 
al. [2009). On the other hand, the estimated envelope mass of 
3.0 Mq and the magnitude of the near-infrared emission (sim- 
ilar to the magnitudes of the UCH n region, which has a lumi- 
nosity of ~ 1000 Lq) are suggestive of MM1 being possibly 
an intermediate-mass YSO. In particular, from the outflow mo- 
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Monge et al. (120081 . 



mentum rate (P) we can estimate the luminosity of the object 
powering the outflow (assuming one single outflow). Using the 
relation from Hatchell et al. (2007) we obtain Lboi — 400 Lq, 
and with the relation from Bontemps et al. (1996) we obtain 
Lboi = 560 Lq, consistent with a driving source of intermedi- 
ate mass. 

We note that concerning the driving source(s) of the out- 
flow(s), we discard 2M00 14256 as a possible candidate because 
its near-infrared excess is much smaller than that of MM1 (see 
Table [7), and lies slightly off from the line joining the high- 
velocity clumps Ba and Ra, or Bb and Rb. The only aspect sug- 
gestive of 2M0014256 being associated with the star-forming 
region is the elongation seen at 3.2 mm, which matches well 
with the position of the 2MASS source (see Fig. QJ). However, 
such an elongation at 3.2 mm could be tracing dust entrained by 
the outflow, as the elongation and the Ba and Ra lobes are very 
well aligned. If 2M0014256 was really associated with the re- 
gion one would expect some compact emission associated with 
it in the 1 mm continuum maps (which have higher angular reso- 
lution than the 3.2 mm continuum maps), but as this is not found, 
it suggests that the elongation is made of extended emission. 

Another aspect favoring the scenario that MM1 could be the 
driving source of the outflow is that MM1 is associated with cen- 
timeter emission, typically found toward Class 0/1 sources driv- 
ing outflows. Since we only detect the source in the centimeter 
range at 3.6 cm we can only estimate a range of possible spectral 
indices, -0.3 < <?mmi < 1-1, which is consistent with thermal 
free-free emission from ionized gas (the spectral index was cal- 
culated from the images at 6, 3.6, and 1.3 cm made with the same 
Mv-range). The fact that the centimeter source is elongated in a 
position angle of ~ 120 + 30°, almost perpendicular to the CO 
outflow, could be due to contribution from the different subcon- 
densations associated with MM1. Additionally, from the work 
of Shirley et al. d2007l l, we estimated the predicted centimeter 
flux which can be accounted for from the outflow momentum 
rate, and obtained a flux density at 3.6 cm of ~ 0.3 mJy. Since 
we measured a flux density at 3.6 cm of 0.17 mJy, we conclude 
that all the centimeter emission associated with MM1 can be ac- 
counted for through ionization by shocks from the outflow. 
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The interpretation of the high-velocity CO bipolar outflow 
and the identification of the driving source(s) are not straight- 
forward. A possibility is that clumps Ba and Ra are tracing an 
outflow with a position angle of 45° powered by one of the 
subcondensations of MM1, maybe MM1-S, as it falls exactly 
at the center of symmetry of clumps Ba and Ra (although it is 
less massive than MMl-main); and that Bb and Rb are tracing 
a second outflow with a position angle of 20°, possibly pow- 
ered by MMl-main, as it shows an elongated structure in the 
east-west direction, almost perpendicular to the Bb-Rb outflow. 
Another possibility is that all the CO emission comes from a sin- 
gle wide-angle outflow, which would be driven by MMl-main 
(since it shows an east-west elongation), and would be excavat- 
ing a cavity, seen at moderate velocities. If this was true, the 
eastern wall of the cavity, corresponding to clump Ba, would 
have very high-velocity emission (up to ~ 35 km s _1 with re- 
spect to the systemic velocity). A similar case in the literature is 
found for G240.3 1+0.07 (Qiu et al. |2009l a massive YSO driv- 
ing a wide-angle outflow, and with velocity in the walls of up to 
27 km s _1 with respect to the systemic velocity. Such a high ve- 
locity in the cavity wall could be explained through precession 
of the outflow axis and episodic mass ejection. Arce & Goodman 
(2001) study the position-velocity relation for episodic outflows 
and show that clumps at different positions from the driving 
source are expected. Each one of these clumps shows a wide 
range of velocities (see e. g., Fig. 2 in Arce & Goodman 1200 II) . 
This can also be seen in the pv-plot of IRAS 001 17+6412 at 20° 
(Fig.fTOl-bottom). In the figure, both the blueshifted and the red- 
shifted sides show a pair of clumps, each one spanning a range of 
~ 20 km s . Thus, both possibilities seem plausible for this tar- 
get, both mutiple outflows driven by different subcondensations 
of MM1, and a single wide-angle episodic outflow. 



5.3. MM2: an intriguing protostellar object 

In Sect. 14.31 we modeled the N2I-F emission toward MM2 as 
a disk structure undergoing infall and rotation. The infall and 
rotation velocities (at 5") can be compared with those derived 
from other protostellar cores. Beltran et al. (2005) model, for the 
high-mass case for example, two massive cores with infall ve- 
locities of 1-2 km s _1 , and rotation velocities of ~ 2 km s _1 , 
at spatial scales of 5000-10000 AU. On the other hand, infall 
and rotation velocities for low/intermediate-mass protostars are 
around 0.1-0.4 and 0.2 km s _1 , respectively, at spatial scales of 
~ 3000 AU ( Beltran et al. [20041 Ward-Thompson et al. [20071 
Carolan et al. 2008). From the infall and rotation velocities at 
the reference radius derived in Sect. l4.3l for MM2, and using the 
assumed infall and rotation power laws (Table [5j, we found an 
infall and rotation velocity of 0.29 and 0.2 km s _1 respectively 
at 3000 AU, which is more similar to the velocities derived for 
low/intermediate-mass cores than to the velocities derived for 
high-mass cores. In addition, the flattened structure was modeled 
with an inner gap of ~ 1". We note that this inner gap must not 
necessarily be a physically real gap. First, it could be an opac- 
ity effect. However, this possiblity can be discarded because the 
opacity of the isolated line (the hyperfine used to fit the model in 
Sect. 14. 3t estimated from the fits is ~ 0.05 (see Table [H, which 
is clearly optically thin. Second, a more plausible option is that 
the N2FF gap is reflecting the depletion of N2H + in the center, 
maybe due to freezing-out of N2I-U due to high densities and low 
temperatures in the center, as found in very young (sometimes 
pre-protostellar) and dense (5 X 10 5 cirT 3 ) regions (Belloche & 
Andre 2004; Pagani et al. |20071 >. High angular resolution obser- 



vations of the continuum millimeter emission would help assess 
the nature of this N2H + gap. 

It is worth noting that the properties of MM2 are quite dif- 
ferent from those of MM1. MM2 has no infrared emission and 
is embedded in a dusty compact condensation and in a dense gas 
clump, which seems to be undergoing rotation and infall mo- 
tions. As there is water maser emission associated with the MM2 
peak (see Fig. QJ), indicative of stellar activity, one would clas- 
sify this source as Class 0, judging from its lack of infrared emis- 
sion and its strong and compact millimeter continuum emission. 
However, YSOs in the Class evolutionary stage are typically 
associated with strong and collimated molecular outflows (e. g., 
Bachiller 1996), while we found no hints of (high-velocity) CO 
associated with MM2. Furthermore, water masers associated 
with YSOs are related to accretion or ejection processes (e.g., 
Furuya et al.|2003), and observational studies suggest a connec- 
tion between water masers and the outflow phenomenon (e. g., 
Zhang et al. [2001). Thus, the lack of outflow emission associ- 
ated with MM2 is intriguing. Different possibilities could ex- 
plain this behavior (lack of outflow in an apparent Class source 
with water maser emission). First, the object could be driving 
an outflow very faint in CO, but which could be well traced by 
other molecules such as SiO (studies of high-mass star-forming 
regions show that some molecular outflows are faint or not de- 
tected in CO while they are strong in SiO, and viceversa, Beuther 
et al. 120041 Zapata et al. 2006). However, a preliminary reduc- 
tion of SiO (1-0) data observed with the VLA (Busquet et al., in 
prep.) show no outflow hints. A second possibility would be that 
the outflow has not been created yet (the protostellar wind would 
still be in the phase of sweeping the ambient material out). The 
caveat with this possibility is that observations seem to indicate 
that outflows appear at the very first stages of the protostar for- 
mation (e. g., Bachiller 1996), and the probability of witnessing 
such a short-lived phenomenon is low. Third, the ejection pro- 
cess in this source could be different from the standard paradigm 
of a strong and collimated outflow associated with a Class 
source. This is found for some objects in massive star-forming 
regions (Torrelles et al. 1200 1 1 120031 >, which show spherical ejec- 
tions traced by water masers. Thus, the true evolutionary stage 
of MM2 remains an open question. A study of the water maser 
emission with the highest angular resolution available would be 
of great help to elucidate the process of accretion/ejection in this 
enigmatic object. 

In order to make a rough estimate of the MM2 luminosity 
we used the water maser luminosity reported by Wouterloot et 
al. (119931 1 toward MM2, of 1.7 x 10" 6 Lq, which we corrected 
to our adopted distance. From a correlation between the water 
maser luminosity and the bolometric luminosity of the YSO as- 
sociated with the maser (Furuya et al. 2003 , 2007), we estimated 
a luminosity for MM2 of ~ 600 Lq, suggesting that MM2 is of 
intermediate-mass. 

5.4. Different objects emerging from the same natal cloud 

The results obtained toward IRAS 00117+6412 reveal that 
the dusty cloud harbours about three intermediate-mass YSOs 
showing different properties. One of the sources seems to be an 
intermediate-mass UCH 11 region with a shell-like structure, lo- 
cated at the eastern border of the dusty cloud, and almost de- 
prived from dust and dense gas, with an estimated luminosity of 
~ 1000 Lq, and an estimated age of ~ 1 Myr. Another source, 
MM1, is deeply embedded within the dusty cloud, has about 
400-600 Lq, and an estimated age around < 10 5 yr. Thus, MM1 
presumably formed after the UCH 11 region. Finally, MM2, of 
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about ~ 600 Lq, remains the most enigmatic object in this re- 
gion, as it seems to be deeply embedded in gas and dust and has 
a water maser associated, but no signs of CO outflow activity. 
This could be a new type of object undergoing a special process 
of matter ejection (such as spherical mass ejection). 

In summary, our observations show that the formation of 
stars within the nascent cluster in IRAS 001 17+6412 seems to 
take place in different episodes. In addition, these observations 
show that the similar initial conditions within a cloud (similar 
dust mass, excitation temperature, and dense gas column den- 
sity, as found for MM1 and MM2) can yield objects with very 
different properties. This indicates that these initial conditions 
may not be decisive in determining some of the properties of the 
YSOs forming within the cloud, such as the ejection properties, 
a result already found at much smaller spatial scales by Torrelles 
et al. (I20TJT1 120031 

6. Conclusions 

In this paper we study with high angular resolution the centime- 
ter, and millimeter continuum, and N2H + (l-0), and CO (2-1) 
emission of the intermediate-mass YSOs forming within a dusty 
cloud, with the goal of assessing the role of the initial conditions 
in the star formation process in clusters. Our conclusions can be 
summarized as follows: 

1 . A UCH n region is found at the eastern border of the dusty 
cloud, with a shell-like structure and a flat spectral index, 
-0.03 + 0.08. The estimated age and mass of the underlying 
star is ~ 1 Myr and ~ 6 Mq. 

2. Deeply embedded within the dusty cloud, we have discov- 
ered a millimeter source, MM1, associated with a 2MASS 
infrared source, which is driving a CO (2-1) powerful and 
collimated high-velocity outflow, oriented in the southwest- 
northeast direction. The mass derived from the millimeter 
continuum emission for MM1 is ~ 3 Mq. MM1 is embed- 
ded within a ridge of dense gas as traced by N2H + , which 
seems to be rotating roughly along the outflow axis. MM1 
is associated with centimeter emission, whose spectral index 
is compatible with an ionized wind, and at 1.2 mm splits up 
into different subcomponents when observed with an angu- 
lar resolution of < 1". From the derived outflow momentum 
rate, we estimated a luminosity for MM1 of 400-600 Lq. 
Thus, MM1 seems to be a Class 0/1 intermediate-mass YSO 

3. About ~ 15" to the south of MM1, our observations revealed 
a dust compact condensation, MM2, lying in a dark infrared 
region, associated with water maser emission and a dense 
core traced by N2H + emission. The mass from the dust emis- 
sion is ~ 1.7 Mq, and the N2FF excitation temperature and 
column density are similar to the ones derived for MM1 . The 
dense core in MM2 is rotating in the same sense as the ridge 
associated with MM1 and seems to be undergoing infall. We 
modeled the MM2 dense core as a disk-like structure with an 
inner radius of ~ 1" and an outer radius of ~ 5", with a ro- 
tation velocity in the outer radius of ~ 0.6 + 0.2 km s _1 , and 
an infall velocity at the same radius of ~ 0.17 + 0.02 km s _1 . 
The non-detection of CO at any velocity toward MM2 makes 
this object intriguing. 

4. Although MM1 and MM2 formed within the same cloud and 
have similar dust and dense gas emission, their properties, 
specially concerning the ejection phenomenon, seem to be 
different and could be indicating that the initial conditions in 
a cloud forming a cluster are not the only agent determining 
the properties of the members of the cluster. 
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Appendix A: Emission filtered out by an 
interferometer 

In this appendix we describe the estimation of the fraction of 
flux filtered out by an interferometer. Consider a bidimensional 
Gaussian source with the flux density S v and the half -power di- 
ameter D. Its intensity, I(x, y), can be expressed as 



I(x,y) = 



41n(2) 



7tD 2 



5 v exp 



-41n(2)(x 2 + y 2 ) 



D 2 



(A.l) 



The visibility, V(u, v), of the Gaussian source corresponds to 
the fourier transform of the intensity (Eq. lA.U . 



V(u, v) = S y exp 



-ln(2) (u 2 + v 2 ) 



[(21n(2))/(;rD)] 2 J 
and thus, the half -power (m, v) radius, ri/2, is 
2 ln(2) 



n/2 



nD 



or in practical units, r\p_ in kA and, D in arcsec 

-l 



n/2 



kA 



91.019 



D 



(A.2) 
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From Eq. IA.4I we can estimate the largest structure, 6>las, 
sensitive to an interferometer with the shortest baseline being 



(A.5) 



From Eq. IA.2l and lA.5l we can estimate the fraction of corre- 
lated flux, V(u m i n )/S v , of a source with half -power diameter D, 
when observed with an interferometer with a shortest baseline of 



0LAS 


= 91.019 
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arcsec 
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with u n 



= exp 



- ln(2) 



V91.019 
in kA and D in arcsec. 
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